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Abstract 

We study the dispersion relation of the excitations of a dilute Bose- 
Einstein condensate confined in a periodic optical potential and its 
Bloch oscillations in an accelerated frame. The problem is reduced to 
one-dimensionality through a renormalization of the s-wave scatter- 
ing length and the solution of the Bogolubov - de Gennes equations 
is formulated in terms of the appropriate Wannier functions. Some 
exact properties of a periodic one-dimensional condensate are easily 
demonstrated: (i) the lowest band at positive energy refers to phase 
modulations of the condensate and has a linear dispersion relation 
near the Brillouin zone centre; (ii) the higher bands arise from the su- 
perposition of localized excitations with definite phase relationships; 
and (iii) the wavenumber-dependent current under a constant force 
in the semiclassical transport regime vanishes at the zone boundaries. 
Early results by J. C. Slater [Phys. Rev. 87, 807 (1952)] on a soluble 
problem in electron energy bands are used to specify the conditions 
under which the Wannier functions may be approximated by on-site 
tight-binding orbitals of harmonic-oscillator form. In this approxi- 
mation the connections between the low-lying excitations in a lattice 
and those in a harmonic well are easily visualized. Analytic results 
are obtained in the tight-binding scheme and are illustrated with sim- 
ple numerical calculations for the dispersion relation and semiclassical 
transport in the lowest energy band, at values of the system parame- 
ters which are relevant to experiment. 
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statistical theory: ground state and elementary excitations; Bloch oscilla- 
tions. 
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1 Introduction 



The dynamics of cold atomic vapours trapped in a periodic external po- 
tential, as created by means of a far-dctuncd standing wave of light, has 
become an active area of research. Studies of ultracold atoms in such a 
confinement have revealed a number of behaviours which are well known in 
solid state physics: Wannier-Stark ladders [1, 2], Bloch oscillations [3, 4] and 
Landau- Zener tunneling [5]. More recently, the observation of interference 
between macroscopic quantum states of Bose-Einstein-condensed vapours of 
®^Rb atoms confined in a vertical array of optical traps has been reported [6]. 
Interference between the condensates in different lattice wells manifests itself 
in falling drops, which are emitted with a period determined by the Bloch 
oscillations induced by the gravitational field and are interpreted as coherent 
matter-wave pulses. 

The properties of a Bose-Einstein condensate (BEC) in a regular array 
of optical wells have been studied theoretically and numerically by a num- 
ber of authors [7 - 13]. Especially relevant for our present purposes is the 
work on an interacting BEC in an infinitely extended lattice potential by 
Berg-S0rensen and M0lmer [7], who demonstrated band structure and Bloch 
oscillations by solving numerically the Bogolubov - de Gennes equations for 
the excitations and the time-dependent Gross-Pitaevskii equation for motion 
under a constant applied force. They also gave analytic results for the exci- 
tation spectrum in the limit where the periodic potential can be treated as 
a perturbation. 

From solid state theory a tight-binding (TB) approach offers the possibil- 
ity of an analytic study of the same phenomena in the opposite limit, when 
the lattice potential is so strong compared with all other energy scales as 
to allow the use of on-site orbitals for the approximate construction of the 
Bloch states. In fact, the proper theoretical tools are the so-called Wan- 
nier functions, which are defined through the Fourier transform of the Bloch 
eigenfunctions in the momentum representation. In the particular case of 
one-dimensional (ID) motion, the energy bands are non-degenerate and the 
states in any given one of them are described from a single Wannier function, 
by repeating it at each lattice site with the appropriate phase relationship 
between different sites. The TB scheme may then be realized by replacing 
the Wannier function by an on-site " atomic" orbital, under suitable limiting 
conditions which may be valid for the low-lying energy bands. 
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While the Wannier functions for motion in a generic periodic potential 
are usually not easy to calculate, it turns out that for the specific problem 
of ID periodicity posed by the confinement of a BEC in an array of optical 
wells one can take advantage of the results contained in a classic paper by 
Slater [14] on a soluble problem in energy bands for electrons. Slater's work 
also discusses the conditions under which the Wannier functions for the low- 
lying energy bands can be accurately approximated by on-site orbitals, which 
in this case are the harmonic-oscillator wave functions associated with the 
parabola describing the bottom of each well. 

The plan of the paper is as follows. In section 2 we first model a di- 
lute BEC confined in a periodic optical potential as a ID dynamical system 
via inclusion of the transverse size into a renormalized contact interaction. 
We then formulate the problem of ID band structure in terms of Wannier 
functions, the only conceptual difference from the analogous problem of ID 
motion of electrons in a periodic potential being that for condensed bosons 
there is no gap between the condensate ground state and the excitation spec- 
trum. In the same section we also formulate the problem of Bloch oscillations 
of the condensate in a single band in terms of the Fourier transform of the 
corresponding Wannier function. 

The reduction of the general formalism to the TB approximation is ef- 
fected in section 3, after recalling Slater's conditions for its validity and 
extending them to an interacting BEC. The TB scheme displays in an ex- 
phcit manner how the energy levels of a BEC in a single well are shifted 
and spread into bands by the periodic replica of the well potential. We give 
special attention to (i) the lowest-lying energy band of the BEC, which con- 
tains its wave-like phase excitations, and (ii) the connection that the TB 
approximation exhibits between the excitations in the next two higher bands 
and the dipole ("sloshing") and quadrupole excitations in an effectively ID 
harmonic well. The discussion suggests how excitations in these two bands 
may be observed. We conclude this section by simple numerical illustrations 
of the dispersion relation in the lowest band and of the current induced by 
a constant external force, for values of the system parameters which are rel- 
evant to the experiments of Anderson and Kasevich [6]. Finally, section 4 
ends the paper with some concluding remarks. 
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2 Band structure and Bloch oscillations of a 
BEC from Wannier orbitals 

We consider a dilute vapour of atoms of mass m in a Bose-condensed state at 
zero temperature, which interact via elastic collisions described by an s-wave 
scattering length a and are confined in the optical lattice created by a laser 
field. The lattice can be modelled as a transverse Gaussian potential times 
a periodic potential along the z direction, namely the confining potential is 

Ui{r,z) = [/oexp(-r^/rfj sin^(/cz) . (1) 

Here Uq is the well depth, scaling linearly with the intensity of the laser 
beam, rib is the transverse size of the beam and k = 27r/A its wavenumber. 
The lattice period is c? = A/2. 

In this work we reduce the above problem in cylindrical symmetry to an 
effective ID problem. To this end we freeze the transverse motions of the 
BEC and renormalize the interactions by a method proposed by Jackson et al. 
[15] to treat an elongated BEC, in which the coherence length as introduced 
by Baym and Pethick [16] is much larger than the transverse size. This 
reduction will be inapplicable, however, to high-energy excitations, when 
the excitation energy becomes comparable with the transverse confinement 
energy and transverse motions can therefore be induced. 

Under the conditions indicated just above, the transverse confinement 
modifies the mean-field interactions by a factor depending on the axial den- 
sity ct{z), the factor being / dxdy \G{x,y; (y{z))\* where G is the transverse 
part of the wave function in the xy plane. Using a harmonic approximation 
for the radial part of the optical potential, we obtain in the present problem 
an effective scattering length per unit area given by 



where uq is the number of atoms per lattice well and E^^ = h"^/ {m(P'). Notice 
that the energy differs by a factor eight from the recoil energy Erecou = 
fi^k"^ /2m. A similar renormalization of Uq in Eq. is negligible, since 
rib > A. 

The effective ID Hamiltonian H that we treat thus is 

H = j dzi)\z)A{z)^{z) + ^9 1 dzip\z)ip\z)iP{z)ip{z) (3) 
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where ipi^) is the field operator, g = A-nfi^a/m and A(z) = — (/i^/2m)V^ + 
Up{z) — fi, with /i the chemical potential and Up{z) = Uosin^i^kz). In the 
following all energies are referred to the bottom of the lattice potential. 



2.1 Lattice symmetry and Bose statistics 

Let us summarily go through the standard Bogolubov procedure [17] for 
treating the boson Hamiltonian (|^). We make the Ansatz tp{z) = (f){z)+ilj{z), 
which separates the coherent condensate function (f){z) from the fiuctuating 
part 'ip(z) in the field operator. Upon substituting this Ansatz in Eq. (||) and 
requiring that linear terms in tp{z) vanish, one obtains the Gross-Pitaevskii 
equation for (j){z): 

'A{z)+g<j)\z)]<P{z)=0. (4) 

We are taking (j){z) as real without loss of generality, with normalization to 
unity over the lattice period. 

To quadratic terms the Hamiltonian in the fiuctuation operators is diag- 
onalized by the canonical transformation 

a 

= E[<(^)4-<Wa«] ■ (5) 

a 

In Eq. d^) the suffix a denotes an appropriate set of quantum numbers and 
the functions {ua,Va} satisfy the Bogolubov - de Gennes equations, 

L{z)u^{z) - Hi{z)v^{z) = EaU^{z) , 
L{z)v^{z)-Hj{z)u^{z) = -E^v^iz). (6) 

Here, Hj{z) = g(f?{z), L{z) = A{z) + 2Hj{z) and the E^s are the (real) exci- 
tation energies that we aim to determine. We explicitly point out that with 
this form of Hj{z) the interactions enter the determination of the excited- 
state wave functions only through the density distribution of the condensate 
ground state. The same assumption has been made in dealing with the exci- 
tations of condensates in harmonic traps [18] and has led to good agreement 
with experiment for the eigenfrequencies of the quadrupolar surface modes. 

For a BEG in a ID periodic potential the quantum numbers a are the 
band index n and the lattice quasi-momentum q. From symmetry under time 
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reversal we have that, if {unq, Vnq} is a solution of Eqs. (|^) corresponding 
to energy Enq , then {vnq,Unq} is a solution corresponding to energy —Enq 
[19,20]. We may for convenience call these the particle-like and the hole- like 
solution, respectively, and limit ourselves to calculate spectrum at positive 
energy. As for the homogeneous Bose fluid we expect gaplessness, since in 
its ground state the whole condensate is in a single energy level placed at 
the chemical potential. Therefore, the excitations of lowest energy (referred 
to the chemical potential) are phonons, i.e. they have a linear dispersion 
relation. 

It is convenient to introduce solutions Znq±{z) which have definite sym- 
metry under time reversal [7], by setting Znq±{z) = [unq{z) ± Vnq{z)] / V2. 
From equations (^) we have 

L^{z)Znq+{z) = EnqZnq-{z) (7) 

and 

L+{z)Znq-{z) = EnqZnq+{z) , (8) 

where L±{z) = L{z) ± Hj{z). Hence, 

L^iz)L^{z)Z^,^{z) = El^Z^,^{z) . (9) 

This equation shows that the functions Znq+{z) and Znq^{z) are the right 
and left eigenvectors of the non-hermitean product operator L^{z)L_{z). 

The ground state of the condensate can be taken at g = in the low- 
est band (n = 0, say) and special care is needed in handling the functions 
Zqq±{z), in order to ensure gaplessness. This point is discussed for a periodic 
condensate in Appendix A, following the arguments given in Refs. [19] and 
[20]. Here we only need to note that an evaluation of the functions Zng-{z) 
is not necessary for the determination of the excitation energy eigenvalues. 
These are calculated by solving Eq. @ for Zng+{z), which reads 

L+{z)L4z)Zng+{z) = E^Zr^q+iz) . (10) 

This equation applies also to the state Zoo+ (z) at zero energy (relative to the 
chemical potential). Of course, the function Zqq+{z) is simply proportional to 
the function 0(z) introduced in Eq. (^. The solutions of Eq. (|I0|) must have 
the Bloch symmetry, which imposes definite g-dependent phase relationships 
between different wells. 
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2.2 Expansion in terms of Wannier functions 

It is well known from solid state theory [21] that any Bloch orbital can be 
expressed as a superposition of the so-called Wannier functions centred on 
the lattice sites. We thus write 

Znq±{z) = ^= Yl exp{ilqd)wn±{z - Id) (11) 
ViV I 

where the index / runs over all sites of the lattice. Of course, the Wan- 
nier functions Wn±{z) have the same parity under time reversal as the cor- 
responding Bloch orbitals Znq±{z). We may also recall that the Wannier 
functions form a complete orthonormal set, i.e. they satisfy the relation 
/ dz w^^{z — ld)wn'±{z — I'd) = 6nn'Sii'- Orthogonality between different 
wells implies oscillatory tails, falling off more rapidly as overlap decreases. 

Eq. (|TTp expresses a Fourier-transform relationship between Bloch or- 
bitals and Wannier functions, embodying the Bloch translational symmetry. 
The g-dependent Bloch orbitals emphasize the extended nature of the lattice 
states and the periodicity of their energy eigenvalues in momentum space. 
The representation in terms of Wannier functions emphasizes instead that 
the lattice is a periodic assembly of identical wells. The Bloch symmetry of 
the Fourier series in Eq. ([TT|) is ensured by the fact that the Fourier coeffi- 
cients do not involve a separate dependence on position z and on lattice site 
index /. 

In the case of a BEG in a lattice, a general consequence of the Wannier- 
function description is immediately evident from Eq. ([TT|). Since the ground 
state is at the chemical potential and has been taken to lie at g = in the 
bottom of the lowest band, the excitations in this band are wave-like modu- 
lations of the condensate phase and their energy vanishes linearly in the limit 
g — i> 0. This fact can be traced back to the existence of a Goldstone phonon 
- a collective mode without restoring force [19] resulting from spontaneous 
symmetry breaking (see the discussion in Appendix A). 

Similarly, for the higher excitation bands Eq. (0) shows that each of 
them is a superposition of localized excitations (described by Wn+{z)), which 
is constructed by imposing definite phase relationships between the various 
sites as dictated by the translational lattice symmetry. This theorem is valid, 
of course, within the Bogolubov approach and has been derived through the 
reduction to one-dimensionality. 
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In practice, the calculation of Wannier functions is in general not an easy 
task. In the limit in which the Wannier functions become almost confined 
inside each lattice well (which may be true for low-lying states), they can 
be expanded in a limited number of "atomic" orbitals. We shall discuss in 
the next section the conditions under which one may replace the Wannier 
functions for a BEC in a periodic optical potential by single-well orbitals and 
thereby obtain explicit solutions in the TB limit. 



2.3 Semiclassical transport in band states 

We proceed to express through the Wannier functions the current carried 
by a periodic Bose-condensed system subject to a constant force F in the 
semiclassical regime. We restrict ourselves to Bloch oscillations in the lowest 
band, described by the Wannier function wo+{z). Through Eq. (p!T| ) taken 
at g = 0, this function is immediately related to the ground state (f){z) of the 
condensate. 

The quantity of interest is the expectation value p{t) of the momentum 
operator p = —ihd/dz in the state Zog(j)+(-2), where q{t) = Ft/Ti. We intro- 
duce the Fourier transform /o(fc) of Wq^{z), 

1 

/o(^) = ~7k= / dz wqj,{z) e^-p{-ikz) . (12) 

V ZTT "'-oo 

Hence, 

1 

Zoq+{z) = -y== /o(^) exp[i(g - k)ld + ikz] (13) 

y iTxJS I J -co 

and a straightforward calculation yields 

^ nj2il-oodk k\fo{k)\'^cos[ld{q{t) - k)] 
^^^^^^ EiIZodk\fom'cos[ld{qit)-k)] ^ ^ 

This is an exact expression, which holds independently of the confining pe- 
riodic potential and of the interaction between the particles. 

It is easily seen from Eq. ([I^) that (i) p{t) is a periodic function of time 
with the Bloch period Tb = h/{Fd) and (ii) p{t) vanishes at the Brillouin 
zone edges, i.e. at g = ±7i/d or t = ±Tb/2. 
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3 Tight-binding approximation 



The aim of the TB approximation is to relate the properties of the periodic 
system in a lattice potential to the solution of the single (isolated) well prob- 
lem. The results of the method are quantitative as long as the single-well 
solutions at adjacent sites do not overlap much with each other. Namely, 
the validity of the TB approximation is limited by the condition that the 
band-width be smaller than the height of the lattice potential. 

Slater [14] has studied the Wannier functions of the "Mathieu problem" 
defined by the Schrodinger equation for non-interacting electrons moving in 
the periodic potential Up{z) = Uosin^^kz) of present interest. He showed 
that for the first five low-lying energy bands the Wannier functions are well 
approximated by Hermite polynomials under the condition Uo > 2Er, re- 
laxing to f/o > Er/A if only the lowest energy band is of interest. Slater's 
argument is summarized in Appendix B and extended there to an interact- 
ing BEC. In essence the mean-field interactions may be viewed in the present 
context as effecting a renormalization in the height of the external potential, 
thereby shifting Slater's thresholds by the amount 

A[/o = ^ . (15) 

TT 

If the threshold conditions are satisfied, the Wannier functions Wn+{z) at 
positive energy could be approximately replaced by the eigenfunctions of the 
harmonic approximation to the single-well potential U^p\z) = UQsm^{kz) 
with z G [-(i/2, d/2\ and Uq = Uq + MErI-k. 



3.1 Band structure for a periodic BEC 

More generally, the evaluation of the energy bands for a periodic condensate 
in the TB approximation would require a preliminary solution of the on- 
site problem. To this end, we separate the on-site parts L^f^{z) out of the 
operators L±{z) entering Eq. (|^) by writing L±{z) = L^±\z) + AL±{z). The 
detailed definitions of these operators are given in Appendix C. The on-site 
operators L^±\z) determine the single- well orbitals (pn±{z) corresponding to 
single- well energy eigenvalues e„, according to 
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(16) 



Again, the condensate in a single well lies at the chemical potential /i^ (see 
the definitions of L^±{z) in Appendix C). Conjugation between fiuctuations 
in the particle number and in the phase must also be taken care of, as shown 
in Refs. [19] and [20]. 

We define the overlap integrals X„(/) = J dz (p*^_{z)ip^j^{z — Id) and the 
interaction integrals Jnil) = J dz ip'^_{z)A [L+{z)L^{z)] ^n+{^ ~ (see the 
definitions of the AL±{z) operators given in Appendix C). By exploiting Eqs. 
(|16D one easily finds from Eq. (^ the dispersion relation 



An adiabatic argument regarding the switching-on of the (spatial-symmetry 
conserving) mean-field interactions shows that the product ip*^_{z)ipn+{z') 
is even under space inversion. Both series in Eq. ( |17|) are therefore real 
quantities. 

Eq. (0) displays the two main changes induced by lattice periodicity 
on the single-well energy levels, namely (i) a shift of each level relative to 
the single-well value e^, and (ii) the broadening of the level into a band. In 
particular, for the lowest band the shift corresponds to a renormalization of 
the chemical potential from the single-well value /i^ to the value /i, which has 
to be evaluated self-consistently from the condition £'oo = /i. The resulting 
gaplessness introduces in Eq. ([T7|) for n = the linear, acoustic-phonon-like 
dispersion of the excitation energies in the long wavelength limit g ^ 0. 

We can now see how the excitations in the bands with n = 1 and n = 2 
may be related to single well excitations. We make use of the present re- 
placement of Wannier functions by on-site orbitals, which are approximately 
described by the eigenfunctions of the ID harmonic oscillator. Evidently, this 
allows us to identify the on-site excitation leading to the n = 1 band with 
the dipolar sloshing mode of the condensed particles within the well. The 
construction of Bloch orbitals from these localized sloshing motions simply 
amounts to fixing, for each value of the wavenumber q, their relative phases 
according to the lattice translational symmetry (see the discussion already 
given near the end of section 2.2). Similarly, the excitations in the n = 2 band 
may be viewed as resulting from localized, single-well quadrupolar modes 
with relative phases depending on q as imposed by lattice symmetry. Of 




T,i Jnjl) (ixp{ildq) 
Y.i1n{l) exp{ildq) 



(17) 
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course, the approximate validity of such a visuahzation of the collective low- 
lying modes of a periodic BEC is subject not only to conditions relating to 
effective one- dimensionality in a mean-field approach, but also to the restric- 
tions associated with the TB scheme. 

The above discussion also suggests how these collective band modes may 
become accesible to experimental or computational observation. Clearly, at 
q = the dipolar mode can be excited by a rigid shaking of the lattice and the 
quadrupolar one by a modulation of the width of all the wells. Excitation of 
dipolar and quadrupolar band modes at finite q would require similar probes 
at given wavelengths. 

3.2 Dispersion relation and Bloch oscillations in the 
lowest band 

We conclude this section by presenting some illustrative numerical results 
for the lowest band in the TB approximation. According to Eq. ([TSD , in 
the case of repulsive interactions the inequality Uq > (0.25 -|- dd/TT)ER must 
be satisfied for the appropriate Wannier function to be approximated by the 
Gaussian function 

^o^z) = j^jJ^exp{-zy2a') . (18) 

The width a of this Gaussian is determined in a first approximation by the 
curvature of the bottom of the single well. A more accurate estimate of a, 
which takes account of the interactions between the particles in the BEC, is 
given in Appendix C. 

In fact, for the calculation of the dispersion relation it is more convenient 
to revert back to Eq. (p!0| ) and directly evaluate the integrals 

Kil) = jdz ^,+ {z) [L+{z)L_{z)] ^,^{z - Id) (19) 

and /(/) = J dz ipo+{z)ipQ_^_{z — Id). This procedure has the advantage of 



involving only the known symmetric function ip()^{z) given in Eq. (18). Since 
Jo(0) ^ 0.99, we obtain 

= fi^ + 2K{1) [cos(grf) - 1] (20) 
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Figure 1: An example of the dispersion relation of phase excitations in a 
BEC confined in an optical potential. We plot the lowest-lying energy band 
as a function of the wave number q (in units of l-n jd) within the ID Brillouin 
zone, for values of the system parameters given by a = 110 Oq, ri\, = 80 fim, 
no = 10\ Uq/Er = 0.25 and a/d = 0.27. 



0|) are referred 



where //^ = K{0) + 2K{1). We recall that the energies in Eq. 
to the bottom of the lattice potential. 

We have calculated the integrals entering Eq. ( pOD with the wave function 
in Eq. (0) for Uq/Eji = 0.25, which corresponds to a value oi a/d = 0.27, 
in the case a = 110 Oq, rih = 80 /im and hq = 10^. These values are 
relevant to the experiment of Anderson and Kasevich [6]. The results are 
K{1) ~ -4 • 10"^ ^2 and ~ 0.1 Er, yielding the value c, ~ 0.02 End/h 
for the speed of phonon-like excitations. Figure 1 shows the corresponding 
dispersion curve over the whole Brillouin zone. 

Using Eq. (0) in Eqs. (|l^) and (|l^) we find the following expression for 
the current induced by a constant external force in states of the lowest band. 



hd J2Zil sin(/grf)exp 


- {ld/2af 




<^^l + 2 Ez=i cos(/gci) exp 


- {ld/2af 



(21) 
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Figure 2: Average momentum of the BEC in the lowest-lying energy band, 
in units of h/d, as a function of the wave number inside the ID Brillouin 
zone, for the case a = 110 Oq, rn, = 80 /im, Uq = 10^ and for various values 
of the parameter Uq/Er. From top to bottom on the RHS of the drawing: 
Uo/Er = 0.25, 0.4, 0.8 and 1.2. 

In practice, in the TB limit the leading term in the numerator is dominant 
and the denominator can be replaced by unity. Namely, 



Figure 2 shows the behaviour of p{q) over the Brillouin zone for various 
values of the parameter a — Uq/2Er in the case a — IIOoq, r/b = 80 //m and 
no = 10^. 

4 Concluding remarks 

In summary, and as already discussed in the work of Berg-S0rensen and 
M0lmer [7], the properties which derive from lattice symmetry for the wave 




(22) 
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functions and the energy eigenvalues of electrons moving in a periodic ID 
potential also hold for a Bose-condensed system. However, Bose statistics 
and spontaneous symmetry breaking imply that there is no gap between the 
condensate ground state and the excitation spectrum. This latter property 
can be met in the energy spectra of electrons in semiconductors only as a 
consequence of accidental degeneracies induced by specific space symmetries, 
a well known instance being provided by graphite [22]. For a periodic Bose 
condensate the gaplessness property, together with the quadratic structure 
of the energies as expressed in Eqs. (H) (see section 2.1), yields on general 
grounds the linear dispersion relation Eoq = Cs{\qd\) + 0{q^d?) for the lowest 
energy band. 

We have seen that the specific form of the optical potential acting on a 
BEC from a laser field admits two important simplifications in the band- 
structure problem. Firstly, the problem for the low-lying bands can be 
mapped into a one-dimensional one by accounting for the transverse cross- 
section of the BEC through a renormalization of the s-wave scattering length. 
This mapping has been shown to be useful in relation to the experiments 
performed by Anderson and Kasevich [12]. Secondly, the ID problem away 
from the weak- confinement limit is most naturally formulated in terms of the 
Wannier functions, which have already been studied for the Mathieu prob- 
lem by Slater [14] . The use of Wannier functions, as demonstrated in section 
2.2, provides an exact formulation of the lattice problem and immediately 
shows that (i) the lowest excitations of a condensate in a periodic potential 
are wave-like phase modulations, and (ii) the higher bands arise from the 
superposition of localized excitations with definite phase relationships. 

An analytic study of semiclassical transport in a periodic condensate sub- 
ject to a constant external force becomes straightforward with the use of 
Wannier functions. The expectation value of the momentum operator is pe- 
riodic with the Bloch period and vanishes at the Brillouin zone edges, as 
shown in section 2.3. 

The work of Slater [14] allows one to specify the conditions under which 
the Wannier functions for the low-lying bands of the BEC in an optical pe- 
riodic potential, and in particular that for the lowest band containing its 
phase excitations, can be accurately approximated by harmonic-oscillator 
wave functions. These functions are simply described by the curvature of 
the on-site well as renormalized by the interactions in the BEC and allow an 
explicit, though approximate, connection to be made between band modes 
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and single-well modes. They also allow very simple calculations on a peri- 
odic BEC, as illustrated in section 3 for the dispersion relation in the lowest 
energy band and for the average wavenumber-dependent momentum in cor- 
respondence to a particular set of system parameters which is relevant to 
experiment. 
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Fisica della Materia through the Advanced Research Project on BEC. We 
thank Professor G. ladonisi for bringing Slater's work to our attention. 

A The zero-energy solutions for the lattice 
problem 

The equation obeyed by the wave function Zoq+{z) at g = 0, that is Zqq^^z) oc 
(f){z), is 

L4z)Zoo+{z) =0 . (23) 

We have chosen Eqq = relative to the chemical potential fi of the condensate 
(see the definition of A{z) immediately below Eq. (|^)). 

In view of the macroscopic occupation of this quantum state, one must al- 
low for fluctuations in the number of particles and for the accompanying shifts 
of the phase of the condensate. This is accounted for by adding in the Hamil- 
tonian a "kinetic energy" term aP'^/2, where P = J dzZQQ+{z)[ilj {z)+'iIj'^ (z)] is 
the number operator generating such a "translation" and 1/a plays the role of 
a mass. The canonically conjugate operator Q = i J dzZoQ_{z)[ip{z) —ip'^i^z)] 
is easily seen to be the phase operator. It satisfies the appropriate com- 
mutation relation provided that the weighting function Zqq_{z) obeys the 
equation 

L+iz)ZooA^) = aZoo+{z) , (24) 

with the condition / dzZoo^{z)ZQo+{z) = 1/2. 

From a more formal viewpoint [19, 20], the Bogolubov equations (^) at 
zero energy admit the eigenvector {Zqq^, ^oo+}- A second solution is provided 
by the eigenvector of the squared Bogolubov operator with zero eigenvalue. 
This implies that the result of applying the Bogolubov operator to the second 
solution is proportional to {^oo+j ^oo+} through a constant to be determined, 
leading to Eq. (p^ . In fact, upon differentiation of the Gross-Pitaevskii 
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equation with respect to the total number of particles n one is led again to 
Eq. (|2^), with Zqq_{z) = 2y/nd{^/nZQQj^{z)) / dn and a = 2ndfi/dn. Clearly, 
the constant a is related to the change induced in the chemical potential n 
by a fluctuation in the number of particles. 



B Wannier functions for low-lying bands 

For the ID Mathieu problem (non-interacting particles moving in the poten- 
tial Up{z) = Uosm'^{kz)), Slater [14] shows that the eigenfunctions fn{k) in 
the momentum representation obey the adimensional finite-difference equa- 
tion 

{p' + ls- ev^) Up) - [Up + 2) + Up - 2)] = , (25) 

where p = kd/n, s = SUq/Er and e = ERy/s/Uo. In the low-lying energy 
bands fn{p) is a slowly varying function of its argument, so that Eq. (^ is 
reduced to the differential equation 

-s^^^ + {p'-eV~s)Up) = 0. (26) 

Evidently, this is the Schrodinger equation for a harmonic oscillator, with 
energy levels e„ = 2n + 1 and eigenfunctions /„(k) oc exp(— K^/2)if„(K) 
where k = s~^^^p and Hn{n) are the Hermite polynomials. 

The Wannier functions appropriate to the low- lying energy bands of the 
ID Mathieu problem are, therefore, single-well orbitals with the same func- 
tional form as for the harmonic oscillator. From his study of the exact Wan- 
nier functions of the same problem Slater [14] then showed that these single- 
well orbitals provide an accurate decription of the first five bands under the 
restriction s > 16 i.e. Uq > 2Er. The condition of validity for the first Wan- 
nier function in Eq. (|l^) can be estimated from Slater's numerical results as 

To extend this argument to the Wannier functions of an interacting BEC, 
we view the mean-field interactions as equivalent to a renormalization in the 
height of the external potential, by an amount that we may estimate as the 
average value of g(p'^{z) over a lattice cell. From the normalization condition 
satisfied by (p^z), in the case of repulsive interactions the height of the periodic 
potential is then lowered by an amount g/d. This leads to Eq. (^) in the 
main text. 
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C Details on the on-site problem 



C.l On-site Bogolubov operators 

The operators L± entering the Bogolubov equations (H) and (|) can be spht 
into an on-site part and a residual part as L±{z) = L± (z) + AL±{z), where 

^(s) ^ _^2/2^ ^ ^(s) ^ ^(s) _ ^ ^27) 

= lL'^ + 2H^/^ (28) 

and 

AL_ = AUp + AHj - A/i , (29) 

AL+ = AL_ + 2AHi . (30) 

In Eqs. i^^) Uj,''> is the single-well potential U^'^ = Up{z)e {I - Az^ /£), 

6{x) being the step function, Hj'^^ is the on-site mean-field potential if}**"* = 
gwQ^{z) and fis is the on-site chemical potential, while the A's are differences 
between lattice and single-well values. Namely, 

AUp{z) = Uo sm\kz)e - 1 j (31) 

and 

AHiiz) = g (|Zoo+(^)r - ko+(^)r) • (32) 

C.2 Variational calculation of Gaussian width 

A crucial quantity in the calculations reported in section 3.2 is the width a 
of ipo+{z), since the integrals entering Eq. (^) are determined by the overlap 
in adjacent wells. A first estimate of a can be obtained from the curvature 
of the bottom of the well: this determines a frequency cu = {UoEji/2h'^y^'^ 
such that a = {h/muoY^'^ . A more refined variational estimate can be based 
on a method first proposed by Baym and Pethick [16], including the effect of 
interactions in the condensate. 

We construct the energy functional 

^t*^] = ^1"^' ^ + /^^ W^WI</>WI' + \9jdz |0(^)r (33) 
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and require that it be minimized by (po^{z) when we set W{z) = muP'z^ 12. 
This yields the energy function 

E{a) = + ]mu;'a' + , (34) 

which is minimized when a obeys the condition 

^4 ^ ( A)2 + ^ . (35) 

Evidently, repulsive interactions broaden the Wannier function. A similar 
estimate of a for the n — \ band states shows that the role of the interactions 
is decreasing with increasing energy. 
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